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Triple differential dijet cross-sections in e^p interactions measured with the HI 
and ZEUS detectors at HERA are presented. The data are compared to Monte 
Carlo simulations which differ in their assumptions about photon structure and 
parton evolution. Effects of the resolved processes of longitudinally polarized vir- 
tual photons at HERA are investigated for the first time. 



1 Dijet Production at HERA 

The production of dijet events at HERA is dominated by processes, in which 
a virtual photon, radiated from the electron, interacts with a parton in the 
proton. In the region of photon virtuality Q 2 ^> Aq^, hard collisions of 
the photons do not necessitate the introduction of the concept of the resolved 
photon (as for the real photon) and the process can in principle be described 
by the direct photon contribution alone, in which the photon interacts as a 
whole with partons from the proton. 

The analyses presented here explore the region Aq CD <C Q 2 < E 2 , where 
different theoretical approaches can be used to take into account higher order 
corrections. 

a) LO direct and resolved interactions based on the DGLAP evolu- 
tion equations and parton showers allows the effects of transversallY-aiid also 
longitudinally polarized resolved photon interactions to be studied Era. Since 
the cross section for longitudinal photons vanishes in the photoproduction 
regime due to gauge invariance, and the concept of a resolved photon breaks 
down for Q 2 > E 2 , the only evidence for this phenomena can be observed 
in the region < Q 2 <C Ef. The main difference between the longitudinal 
(7U and transverse (7^) virtual photon arises from the dependence of the 
respective fluxes on y. While for y — > 0, both transverse and longitudinal 
fluxes are approximately same, the longitudinal flux vanishes for y — ► 1. Also 
the dependence of the point-like (i.e. perturbatively calculable) part of the 
photon Parton Distribution Functions (PDF) on Q 2 and E 2 differs - while the 
7y PDF is proportional to \n(Ef /Q 2 ), the 7^ PDF depends on the scale in a 
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typically hadron-like manner. 

b) kf unordered initial QCD cascades accompanying the hard pro- 
cess, employed for example in BFKL or CCFM evolution, can lead to final 
states in which the partons with the largest kt may come from the cascade, 
and not, as in DGLAP evolution, from the hard subprocess. Such events 
may have a similar topology as is observed for the resolved interactions in- 
troduced in the previous, paragraph. This possibility is investigated using the 
CASCADE generator Ha based on the CCFM evolution equation. 

c) NLO calculations with and without the concept of the resolved vir- 
tual photon offer another natural possibility to be examined. Such compar- 
isons are envisaged. 

All of these three approaches include higher order corrections to the LO 
QCD matrix elements; however, each of the approaches treats them differently. 

2 Measurement of Dijet Cross-Section 

The selection criteria defining the HI and ZEUS dijet samples are summarized 
in Table pi: Q 2 denotes the photon virtuality, y is the inelasticity, s is the total 



HI (16.3 pb- 1 , 1999) 


ZEUS (38.2 pb" 1 , 1996-97) 


y/s = 318 GeV 


y/s = 300 GeV 


2 GeV^ < Q 2 < 80 GeV^ 


0.1 < Q 2 < 10 000 GeV^ 


0.1 < y < 0.85 


0.2 < y < 0.55 


gjetWA > g Ge y 

E t > 6 GeV 


Ef 1 1 > 7.5 GeV 
E ] t et2 > 6.5 GeV 


-2.5 < rf et1 ' 2 < 


-3 < rf et1 ' 2 < 



Table 1. Selection criteria of the dijet samples. 



electron-proton centre-of-mass energy, E 1 ^ 1 ' 2 and rf ct1 ' 2 are the transverse 
energy and pseudorapidity of the jet with the highest or second highest E t , 
and ~E t is defined as {E\ ct 1 + Ef t2 )/2. 

The measured data are corrected for detector effects using a bin-to-bin 
method (ZEUS) or Bayesian unfolding (HI). The largest source of systematic 
errors arises from the main calorimeter calibration uncertainty and, in the 
case of HI, also from a model dependence of the detector correction. The 
ZEUS measurement was presented in more detail at the EPS 2001 
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3 Results and Discussion 



The corrected triple-differential dijet cross-section measured by ZEUS as a 
function of Q 2 , E t and x 1 is shown in Fig. 0. A prediction of HERWIG with 
the SaSID parameterization of the 7^. PDF, as well as the direct contribution 
of HERWIG is compared to the data. Since the overall normalization of 
the LO Monte Carlo simulation is to some extent uncertain, the HERWIG 
prediction has been normalized to the highest £ 7 bin (x 7 > 0.75) in the data. 

The normalization is done separately for each (Q 2 , E t ) bin. 

In the region where Q 2 > E t , the data are well described by the direct 
HERWIG component only. Resolved interactions are needed if Q < E t . 
However, even with the 7^ resolved processes included, HERWIG tends to 
underestimate the data in the lowest Q 2 region for x 1 < 0.75. 

This fact is also demonstrated in Fig. || by the ratio of a(x 1 < 
0.75) ja{x 1 > 0.75). The slope of this distribution can be interpreted as a 
suppression of the virtual photon structure with increasing photon virtuality. 

The corrected triple-differential dijet cross-section measured at HI as a 

2 1 1 

function of Q 2 , E t and x 7 is shown in Fig. pL The HI measurement is per- 
formed in a different phase space (see Table ]l]) and Monte Carlo predictions 
are not normalized to the data. 

A comparison of the HI measurement with HERWIG leads to similar 
conclusions as drawn above for ZEUS. In addition, we can see that for the 
highest Q 2 range (25 < Q 2 < 80 GeV 2 ) and x 7 < 0.75, the HERWIG direct 
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contribution almost describes the data in the lowest E t bin, but is too low 
2 

in the highest E t bin. This indicates an importance of the resolved processes 

2 

even at high Q 2 , once the hard scale, E t , is high enough. 

Standard HERWIG with direct and 7^ resolved contributions underesti- 
mates the data. The description is improved by adding 7^ resolved photon 
interactions, which is done using a slightly modified HERWIG with the cor- 
rect longitudinal photon flux and a recent 7^ PDF parameterization 13. As 
demonstrated in Fig. |3| the 7^ resolved contribution is significant, and brings 
HERWIG closer to the measurement. 

On the other hand, a simple enhancement of the PDF of the 7^, in the 
resolved contribution could lead to a similar prediction as the introduction of 
7£ . To eliminate this ambiguity, the dijet cross-section has also been studied 
as a function of Q 2 , x 7 and y, which is shown in Fig. [|. HERWIG is below 
the data, even if the resolved 7^ is added. This may be due to the uncertainty 
of the overall normalization of the LO Monte Carlo prediction. In the region 
where x 1 < 0.75, the slope of the HERWIG prediction depends significantly 
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on whether 72 processes are included or not. The 7^ contributes significantly 
at low y, while it becomes very small compared to 7J at high y. Unlike a 
pure enhancement of 7^ resolved processes by a constant factor, the addition 
of 72 brings the y dependence of HERWIG much closer to the measurement. 

As motivated in Section [l], the measured cross-sections in Fig. || are also 
compared to a prediction of the CASCADE MC program based on the CCFM 
evolution approach. This theoretical concept does not involve any informa- 
tion about the virtual photon structure and involves many fewer free pa- 
rameters for tuning than the usual DGLAP-like MC programs. Neverthe- 
less, CASCADE describes the data well, except for the Q 2 dependence. The 
Q 2 behavior, however, is related to the parameterization of the unintegrated 
PDFs used in the program, which are not yet constrained unambiguously. 



4 Conclusions 

The dijet cross-sections measured as a function of Q 2 , E t , x 1 and Q 2 , x 1 , 

y at HI and ZEUS have been presented. 

In the DGLAP evolution scheme, the importance of the 7^ resolved pho- 

2 

ton interactions is clearly demonstrated in the region where E t > Q 2 . Addi- 
tional 72 resolved photon contributions further improve the agreement of the 
HERWIG 5.9 prediction with the measurement. 

Exploring the CCFM approach, the MC program CASCADE 1.0 gives a 
qualitative description of the measured differential cross-sections; however, the 
Q 2 dependence is not reproduced. On the other hand, the x 1 dependence in 
CASCADE is comparable to the sum of the direct and resolved contributions 
in DGLAP-like MC programs. 
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Figure 1. Triple differential dijet cross-section d 3 a ep /dQ 2 dE t dx-y for the ZEUS data is 
depicted by points, the shaded bands display an uncertainty estimate arising from the 
main calorimeter calibration. The dashed histograms represent the direct part of HERWIG 
only, the full line stands for the sum of the direct and transverse resolved contributions of 
HERWIG. 
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Figure 2. The ratio of cross-sections R = cr(xj < 0.75)/<r(a; 7 > 0.75) from Fig. 
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Figure 3. Triple differential dijet cross-section d 3 a ep /dQ 2 dE t dx-y for the HI data depicted 
by points is compared to predictions of the HERWIG and CASCADE MC programs. 
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Figure 4. Triple differential dijet cross-section d 3 a ep /dQ 2 dxjdy. 
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